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Abstract 

The discrepancy between the measured and theoretical production cross section of b quarks at 
the Tevatron can probably be explained by the recently proposed scenario of light gluinos of mass 
12 — 16 GeV and light sbottoms of mass 2 — 5.5 GeV. In this scenario, we study a related process 
at the Z pole, Z — > bb\g + bb\g followed by g — ► bb\ / bb\. The hadronic branching ratio for 
this channel is (1 — 3) x 10~ 3 , which is of order of the size of the uncertainty in Rt,. We find that 
a typical event consists of an energetic prompt bottom-jet back-to-back with a "fat" bottom-jet, 
which consists of a bottom quark and two bottom squarks. Such events with a 10 -3 branching ratio 
may affect the measurement of R^; even more interesting if the "fat" bottom jet can be identified. 
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I. INTRODUCTION 



There has been a persisting discrepancy that the measured cross section of hadronic 
production of b quarks measured by both CDF and D0 collaborations is about a factor 
of two larger than the prediction in perturbative QCD with the most optimal choice of 
parameters, such as 6-quark mass (mj) and the factorization scale /z, tuned to maximize the 
calculated rate. 1 Recently, Berger et al. || interpreted the discrepancy in the scenario 
of light gluinos and light sbottoms. Light gluinos of mass between 12 — 16 GeV are pair- 
produced by QCD qq and gg fusion processes, followed by subsequent decays of gluinos, 
g — > bb* I bbi, where the sbottom has a mass 2 — 5.5 GeV. Therefore, in the final state there 
are bb + bib*, and the sbottoms either remain stable or decay into other light hadrons (e.g. 
via i?-parity violating couplings) and go into the 6-jets. Gluino-pair production thus gives 
rise to inclusive 6-quark cross section. The mass range of gluino is rrig = 12 — 16 GeV and 
sbottom mr = 2 — 5.5 GeV. Such masses are chosen so that both the total cross section 
and the transverse momentum spectrum of the 6-quark are reproduced. Before Berger et 
al.'s work, there have been some studies in the light sbottom and/or light gluino scenario 
|J. However, such a scenario cannot be ruled out, unless there exists a sneutrino of at most 
1-2 GeV. 

Such a scenario easily contradicts other experiments, especially, the Z°-pole data because 
of the light sbottom. However, it can avoid the Z-pole constraints by tuning the coupling of 
Zbib\ to zero by choosing a specific mixing angle 9b of bi and &r: sin 2 8b = § sin 2 9w, where 
9w is the Weinberg mixing angle. In spite of this, subsequent studies |7|, ^ |9| showed that 
such light gluino and sbottom will still contribute significantly to Rb via one-loop gluino- 
sbottom diagrams. In order to suppress such contributions, the second b 2 has to be lighter 
than about 180 GeV (at 5a level) with the corresponding mixing angle in order to cancel 
the contribution of bi in the gluino-sbottom loop contributions to Rb- Although Berger et 
al.'s scenario is not ruled out, it certainly needs a lot of fine tuning in the model. In other 
words, instead of saying this scenario is fine-tuned, we can say that so far the light gluino 

1 Refs. [|[ H argued that if the most up-to-date B fragmentation function is used the observed excess can 
be reduced to an acceptable level. Fields 0] interestingly pointed out that correlations between the b and b 
can be used to isolate various sources of production, especially, in his study he included the fragmentation 
of gluon and light quarks. 
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and light sbottom scenario is not ruled out. It definitely deserves more studies, no matter 
whether it was used to explain the excess in hadronic bottom-quark production or not. 

The light gluino and light sbottom scenario will possibly give rise to other interesting 
signatures, e.g., decay of Xb info the light sbottom JTU], enhancement of ttbb production 
at hadron colliders |TTJ, decay of T into a pair of light sbottoms [ |T2[] , and affecting the 



Higgs decay [|T3J. In a previous work ]TJ], we calculated the associated production of a 



gluino-pair with a qq pair and compared to the standard model (SM) prediction of qqbb at 
both LEPI and LEPII (here q refers to the sum over u,d,c,s,b). We found that at LEPII 
the qqgg production cross section is about 40 — 20 % of the SM production of qqbb, which 
may be large enough to produce an observable excess in qqbb events [14]. This is rather 
model- independent, independent of the mixing angle in the sbottom, and is a QCD process. 

In this work, we present another interesting channel in Z decay in the light gluino and 
light sbottom scenario: 

Z -> bb\g + bhg ; followed by g bbl/bh . (1) 

Since the gluino is a Majorana particle, so it can decay either into bb\ or bb\. The final 
state can be bbb*b*, bbbibi, or bbbibl. This channel, unlike the one mentioned in the previous 
paragraph, depends on the mixing angle of bi and &r in the bb\g coupling. 

The hadronic branching ratio of this channel will be shown to be (3.4 — 2.5) x 10~ 3 
for sin20 6 > and (1.4 - 1.1) x 10" 3 for sin 26 b < 0, and for m § = 12 - 16 GeV and 
mi = 3 GeV, which is of order of the size of the uncertainty in R b . The process is the 
super symmetric analog of Z — > bbg, but kinematically they are very different because of 
the finite mass of the gluino and sbottom. A typical event consists of an energetic prompt 
bottom-jet back-to-back with a "fat" bottom-jet, which consists of a bottom quark and two 
bottom squarks. If such events cannot be distinguished from the prompt bb events, they 
may increase the R b measurement {Rl* p = 0.21646±0. 00065 fL5|) with a hadronic branching 



ratio of (1—3) x 10 -3 . If the "fat" bottom jet can be distinguished from the ordinary bottom 
jet, then this kind of events would be very interesting on their own. It is a verification of 
the light gluino and light sbottom scenario. Furthermore, if the flavor of the bottom quarks 
can be identified, the ratio of bb :bb : bb events can be tested (theoretically it is 1 : 1 : 2) ||. 

The paper is organized as follows. In the next section, we present the calculation, in- 
cluding the decay of the gluino into bb\ or bb\. In Sec. Ill, we show the results and various 
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distributions that verify the "fat" bottom jet. We conclude in Sec. IV. There is an analog 
in hadronic collisions, pp — > bb\g followed by g — > bb\ / bb\. Thus, it also gives rise to two 
hadronic bottom jets. However, in hadronic environment it is very difficult to identify the 
"fat" bottom jet. We believe it only gives a small correction to the inclusive bottom cross 
section. 

II. FORMALISM 

The interaction Lagrangian among the bottom quark, sbottom, and gluino is given by 

£ D ^g s [b\J«(sm9 b P L + cos9 b P R )T^ + h.c. ] , (2) 

where the lighter sbottom b\ is a superposition bi = sxn9 b bi + cos 9 b b R of the left- and 
right-handed states via the mixing angle 9 b . As mentioned above, the vanishing of the Zb\b\ 
coupling requires Ql sin 2 9 b + g R cos 2 9 b = 0, where 9l = —\ + \ sin 2 9w and g R — | sin 2 0w- 
It implies sin 2 9 b = | sin 2 9w- 

A. Primary Production 

Even a perfect cancellation in the amplitude Z — > b\b\, the Z boson can still decay at 
tree level into bb\g (or its conjugated channel) as shown in Fig. [I]. The Feynman amplitude 




FIG. 1: The Feynman diagram for the process Z — > bb\g. 

is 

M = V2g 8 g z u(b) jt z (g L P L + g R P R ) - / + ™ b (sin 9 b P R + cos 9 b P L ) Tg v(g) , (3) 

p z — m b J 

where Pl, R — (1 T 7 5 )/2, gz = 9ij cos^vi/, and a correspond to the color indices of 

the final-state particles b, b\ and g, respectively. We can tabulate the complete formula of 
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the transitional probability, summing over the initial- and final-state spin polarizations or 
helicities, and colors, as 

16o 2 a 2 

E l-^l 2 = ( 2 + m.m-gN, sm28 b + m 2 b N 2 ) , (4) 

\P ~ m b) 

No = (gl sin 2 9 b + g\ cos 2 9 b )[Ag ■ p b ■ p + p 2 g ■ b (p 2 - m 2 - 2s) / 's + 2g ■ p p ■ Zmf/s] 
N! = 3(p 2 + m 2 )g L g R + {g\ + g 2 R )(p ■ b + 2p ■ Z b ■ Z/s) (5) 
N 2 = Qg L g R g-p+{g 2 R sin 2 b + g\ cos 2 d b ){g ■ b + 2g ■ Z b ■ Z/s) , 
where s = M§. Here the momenta of the particles are denoted by their corresponding 
symbols. We use p to denote the momentum of the virtual b, which turns into g and b\ (i.e. 

p = g + k)- 

One can integrate the exact 3-body phase space to find the decay rate, 

d T(Z^bhtg) = lY,\M\^ d -^§i. (6) 

The scaling variables of the 3-body phase space are defined by 

x b = 2E b /M z , X' b = 2E- K /M Z , and x- g = 2E~ g /M z , (7) 

with the energies Ei measured in the Z rest frame, and x b + x b + Xg = 2. The ratios of the 
mass-squared are 

fi b = m 2 h l M\ , H = m\jM 2 z , and fi 9 = m 2 /M 2 z . (8) 
The region of the phase space is limited by 

1\J~H < %b < 1 + H - H ~ - 1\JWg > ( 9 ) 

X~ b <|(1 - X b + yUb)~ 1 [(2 - X b )(l + /I b + m~ fig - X b ) ± (X 2 b - 4/i 6 )5A5(l + fl b - X b , fl~ b , fig)) 

with the function A(a, b, c) = a 2 + b 2 + c 2 — 2ab — 26c — 2ca. The scalar dot-products can be 
expressed in terms of the scaling variables as 

p 2 = s(l + fi b - x b ) , g ■ b = |s(l - x- b + m - fig - fi b ) 

b-p= ±s(x b - 2fi b ) , g -p= \s{\ - x b - fi~ b + fi § + fi b ) 

The calculation for the charge-conjugated process Z — > bb\g can be repeated in a straight- 
forward manner. Eqs. (0) and (Q) remain valid if we make the substitutions b b, 6* b\. 
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B. Decay of gluino 



Since the gluino so produced will decay promptly into bb* or bb\, the event ends up 
with the final states bbb\b\, bbbib*, or bbbibi. In the minimal hypothesis that the sbottom 
hadronizes completely in the detector, it behaves like a hadronic jet. The final configuration 
includes bb + 2j, bb + 2j, and bb + 2j at the parton level. We will show below that the 
2j most of the time goes together with the softer b, and therefore makes the b look "fat" . 
The complete jet structure requires the full helicity calculation following the decay chain 
Z — > bb\g and g — > bb\ or bb\. Based on Feynman rules for the Majorana fermions, we 
replace v (g) in the above Eq. (^) by 

Ch.l : v(g) — -V2g s T a - - \ + m ^ (sin0 6 P L + cos 6 b P R ) v{b) , (10) 

g z - mj + i\-~grrig 

for the process g —*bb\. Similarly, we replace v(g) in Eq. (||) by 

Ch.2: v(g) — > V2g s T a — ~ { + ™f (sm6 b P R + cos9 b P L ) v(b) , (11) 

g 2 - mj + «r 5 mg 

for the process g —> bb*. We use the narrow-width approximation to calculate the on-shell 
gluino propagator 

1 " S(f-ml), (12) 



(p-rrqy + Tjml m g T g vy 
where g — b + 6* or b + b\. Assuming the gluino only decays into bb\ and bb\, we find the 
decay width of the gluino is 

T g = \(a s /m s ) \*(l,ml/m~,mljrn~) (m? + ml - +2m g m b sm29 b ) . (13) 

Since we have already assumed CP invariance in Eq. @, the event distributions of a pair 
of CP-conjugated variables are the same. 

III. RESULTS 

We first list the input parameters in our study 



2 2 2 2 

m b = 4.5GeV, m^=3GeV, sin# fe = W-sin 6 W , cos6 b = ±W 1 — - sin ^ 
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The scale Q that we used in the running strong coupling constant is evaluated at a s (Q = 
M z /2). 2 

We show in Fig. |^ the partial width of the channel Z — > bb^g + bb\g versus the gluino mass 
irtg for two different sign choices sin26 l {, <0. Numerically, the effect of is not negligible 
at -y/i = Mz- Given that the total hadronic width of the Z boson is 1.745 GeV [|15 



the hadronic branching fraction of the process Z — > bb\g + bb\g is (3.4 — 2.5) x 10~ 3 for 
sin26> 6 > and (1.4-1.1) x 1(T 3 for sin26> 6 < 0, and m- g = 12-16 GeV. Thus, this hadronic 
branching ratio is at the level of, or even larger than, the uncertainty in the Rb measurement 
(i?b Xp = 0.21646 ± 0.00065). If it cannot be distinguished from the prompt bb events, it will 
affect the precision measurement on the bb yield at LEP I. 

In the following, we study the event topology to examine the difference from the prompt 
bb production, which essentially consists of two back-to-back clean bottom jets with energy 
equal to Mz/2. In Fig. ^|, we show the energy distributions, in terms of dimensionless 
variables Xb, x^, Xg, of the prompt b, sbottom, and gluino, respectively. The prompt b has a 
fast and sharp energy distribution as expected, but the gluino and the sbottom have slower 
and flatter energy spectra. We also note that the spectra are different between sin 28b > 
and < 0. These features are very different from the prompt bb production including QCD 
correction, in which both b and b are very energetic and the gluon is quite soft. 

In Fig. [|, we show the energy spectra for the decay products, &dec and &dec, of the gluino. 
Since gluino is a Majorana particle, it decays into either bb{ or bb\. Although there are 
some differences between these two decay modes because of the difference in the coupling, 
in both modes the &dec and fedec are rather soft. We also note that the spectra are different 
between sin 26*6 > and < 0. Therefore, just by looking at the prompt b and the secondary 
&dcc the energy spectra are very different from the prompt bb production. However, if the 
first and the second sbottoms go very close with the secondary 6 d ec and cannot be separated 
experimentally, and the sbottoms deposit all their energies in the detector, then the event 
will mimic the prompt bb event. Thus, it is important to look at the angular separation 



2 The difference in a s between including and not including the light gluino and sbottom in the running 
of a s from Q = Mz to Alz/2 is only 3%. Thus, we neglect the effect of light gluino and sbottom in 
the running of a s . Refs. || [l^] also estimated the effect of including the light gluino in the running of 
a s in their studies. A recent work fill studied the running of a s from low-energy scales such as m T to 
Mz including a light gluino and a light sbottom. However, it cannot rule out the existence of such light 
particles from current data. 
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among the final-state particles. 

We show the cosine of the angles between the primary b and the bdec, between 6d ec and 
&dec> and between b^ ec and b\ in Fig. |5|. Here we only show the spectra for the case sin 29b > 
and gluino decay Ch. 1, because for sin 29 1 > or < 0, Ch. 1 or Ch. 2, the spectra are very 
similar. We can immediately see that the primary b is back-to-back with the secondary b& ec 
from gluino decay. The &dec and 6dec are very much close to each other that the cosine of 
the angle between them is peaked at 0.8 — 0.9. The cosine of the angle between b& ec and 
b\ has a broader distribution, but still peaks in the cos# = 1 region. Thus, we have the 
following picture. The decay products, &dec and &dec, and the primary b\ combine to form a 
wide or "fat" bottom-like jet. This "fat" bottom jet is back-to-back to the primary bottom 
jet, which has an energy close to Mz/2. 

Here we comment on the possibility that the channel that we consider here may affect 
the Rf, measurement, based on two criteria. First, one of the bottom jet in the channel 
under consideration is "fat". If the two sbottoms cannot be separated from the bottom, 
the resulting bottom jet will just look like a fat bottom jet and may affect Rb. Second, 
whether the energy in this fat bottom jet equals to half of the Z mass or not. As mentioned 
by Berger et al. |J, the sbottom can either decay into light hadrons or escape unnoticed 
from the detector. If the sbottoms escape the detection (which means that they do not 
deposit enough kinetic energy in the detector material for detection), the fat bottom jet 
would have an energy much less than Mz/2. The final-state would be two bottom jets (one 
energetic and one much less energetic) plus missing energy, and thus would not affect Rb- 
Neverthesless, this is a very interesting signal on its own. On the other hand, if the sbottoms 
deposit all their kinetic energy in the detector, the measured energy of the fat bottom jet 
would be close to Mz/2. In this case, it may affect the measurement of Rb- In fact, it would 
increase Rb- But if the fat bottom jet could be distinguished from the normal bottom jet, 
the present channel is also interesting on its own. According to a study on the light gluino 
17f , a sbottom of mass 2 — 5.5 GeV, if similar to gluino, will likely deposit most of its kinetic 



energy in the detector. If this is the case the signal would be two back-to-back bottom jets, 
one of which is "fat" or wide, with no or little missing energy. 



S 



IV. CONCLUSIONS 



We show that the light-sbottom-gluino scenario predicts the production of 666*6*, 66616*, 
and bbbibi at the Z pole, with a branching fraction of order of 10~ 3 , depending on the 
gluino mass and the sign of the mixing angle. The event topology is very different from the 
prompt 66 production. Depending on whether the sbottoms deposit little or almost all of 
their energies in the detector, the signal would be very different. If the sbottoms escape the 
detector unnoticed, the final-state would be two bottom jets (one energetic and one much 
less energetic) plus missing energy. On the other hand, if the sbottoms deposit all their 
kinetic energy in the detector, the final state will be two bottom jets, one of which is fat. 
In this case, it may increase the measurement of Rf,. But if the fat bottom jet could be 
separated from the normal bottom jet, it is a distinct signal. These two kinds of signals may 
well be hidden in the LEP I data, waiting for deliberate search. 

One special feature of the Majorana nature of the gluino predicts a ratio of 1:1:2 for the 
rates of 66 : 66 : 66 0. However, one needs to look for the charged modes B + B + or B~B~ 
to avoid effects due to B°-B° oscillation. 
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FIG. 3: Normalized energy spectra of the b (x^), sbottom (x^), and gluino (x g ) in the decay 
Z bb\g. (a) sin26> 6 > and (b) sin20 6 < 0. 
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(b) sin29 b < 
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FIG. 4: Normalized energy spectra Xb dcc and , in which the bottom and sbottom are the 
subsequent decay products of the gluino. (a) sin26^ > and (b) sin 26*6 < 0. Here "ch. 1" and 
"ch. 2" refer to the decay channels of the gluino in Eqs. (|l0| ) and (0), respectively. 
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FIG. 5: Normalized spectra of the cosine of the angles between various pairs of final-state partons 
in the decay process Z — > bb\g, followed by g — ► 661. Here &dec and &dec denote the decay products 
of the gluino. 
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